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Abstract. When a theory shall be described at all scales, it is necessary to start from its elementary degrees of freedom. 
Herein, one possible chain of steps for this purpose will be briefly outlined for the example of a gauge theory, like QCD. 
Starting with the elementary constituents, gluons and quarks, step by step the final observables, physical states, will be built 
up. This process is based on the elementary correlation functions, and uses a combination of both numerical lattice calculations 
and functional continuum methods. While the process uses a fixed gauge at each intermediate step, the final observables are 
gauge-invariant. 
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WHY A GAUGE-FIXED APPROACH? 

In principle, it is possible to calculate any observable 
quantity in a gauge theory without fixing a gauge. How- 
ever, for essentially all practical calculations in non- 
integrable gauge theories it is useful to fix a gauge at 
least for some of the intermediate steps for at least two 
reasons. On the one hand, many calculations simplify 
considerably when choosing a suitable gauge. Especially 
perturbation theory provides a rich set of examples for 
this fact. On the other hand, fixing the gauge permits to 
construct intermediate results, like the correlation func- 
tions of the elementary degrees of freedom. These can 
then be used to construct observables out of them. This 
is, e. g., done in Bethe-Salpeter or Faddeev equations for 
hadronic properties [1]. Furthermore, they can be used to 
exchange results between different methods, e. g., lattice 
and functional methods [2]. The latter case is particularly 
important if a combination of methods is applied to solve 
a certain problem which cannot be solved using a single 
method. A pertinent example is to use functional meth- 
ods to extend lattice results from the narrow energy do- 
main accessible in numerical simulations to the energy 
reach of at least 9 orders of magnitude typical for the 
standard model. 



process of gauge-fixing, due to the so-called Gribov- 
Singer ambiguity. A short description of the inherent 
problems can be found in [1, 3]. An important con- 
sequence of this ambiguity is that in general a pertur- 
batively unique gauge description becomes a family of 
gauges non-perturbatively. 

One such family which has been rather intensively 
studied is the Landau gauge family, and herein in partic- 
ular the subset of so-called first-Gribov-region Landau 
gauges [3]. The reason for the choice of Landau gauge 
is its technical advantages, for details see [2]. However, 
many other gauges have been studied as well, see e. g. 
[1] for an introduction. It should be noted that it is rather 
non-trivial to choose the same gauge with different meth- 
ods beyond perturbation theory, see [2, 3, 4]. In the re- 
mainder, it will be assumed that the construction pro- 
posed in [3, 4] is possible. 

Once a choice of gauge is done, the elementary de- 
grees of freedom are well defined. It is then possible to 
determine their correlation functions, in particular their 
propagators and vertices. This will be done here using 
both lattice and functional continuum methods. 



PROPAGATORS 



GAUGES 

In a gauge theory, the elementary degrees of freedoms 
cannot be described without fixing a gauge. E. g., in 
QCD these are the quarks and gluons. Fixing the gauge is 
merely a choice of coordinate system in the field configu- 
ration space. This choice is often facilitated by the intro- 
duction of additional auxiliary fields, the so-called ghost 
fields. However, gauge theories based on non-Abelian 
gauge groups introduce additional complications in the 



Once the gauge is fixed, it is possible to determine both 
propagators [2, 9] as well as higher vertex functions 
[9, 10], though this can be rather challenging in detail. 
Results for the propagators of the gluon, the ghost, and 
the /3 -function of Yang-Mills theory for various gauges 
and methods are shown in figure 1 . It should be noted that 
the renormalization scheme has been chosen differently 
for the different gauges for convenience, and therefore 
the j3-functions do not agree for the different gauges. 
That this is also possible at the level of vertices is shown 
in figure 2 for the ghost-gluon vertex. 
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FIGURE 2. The single tensor structure of the ghost-gluon 
vertex from lattice calculations in the standard gauge [7] and 
from DSEs in the scaling gauge [8] at the symmetric point, i. 
e., for all momenta of the same magnitude [8, 9]. 



can then be used to obtain physics results. 
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FIGURE 1. The gluon propagator (top panel), the ghost 
dressing function (middle panel) and the jS-function (bottom 
panel, axes rescaled) for lattice calculations using a standard 
gauge [5], its reproduction using Dyson-Schwinger equations 
(DSEs) [2], as well as DSE results for a scaling gauge [2] and 
for a third gauge in lattice calculations, the maxB gauge [4, 6]. 
Concerning the assignment of the gauges, see [3]. 



The results show that it is possible to obtain rather well 
agreeing results with functional and continuum calcula- 
tions, at least on the level of propagators [2]. Thus, it 
is possible to determine the correlation functions, which 



The main aim of this combination of methods is to ob- 
tain results which cannot be easily determined with ei- 
ther method alone. Particular examples are chiral limit 
physics [1 1] or the finite density region in the QCD phase 
diagram [12]. Useful quantities in this context are meson 
masses [1, 11] and order parameters, like the Polyakov 
loop [13, 14, 15]. 

As an example, consider the mass of a meson [11, 
16, 17]. For simplicity take the quenched case. Then, 
the mass of the meson can be determined using the 
Bethe-Salpeter equation [1], which requires further input 
in form of the quark propagator and the quark-gluon 
vertex. This propagator can be determined again using 
the same methods, and will depend on the chosen gauge 
[11, 16, 18]. Results for different gauges for its two 
independent dressing functions are shown in figure 3. 

The second ingredient necessary is the quark-gluon 
vertex. It can be either calculated self-consistently [10], 
or it can be parametrized using the gauge-invariance 
of meson masses to fix the parameters [16]. Fixing the 
vertex using the pion properties with the latter method, 
it is possible to calculate, e. g., the p -meson properties 
[16]. Since the gauge dependence is, up to truncation 
artifacts, canceled between the different contributions 
from the gluon, the quark, and the quark-gluon vertex, 
the p -meson properties then also agree with the physical 
values for physical quark masses [16]. On the other hand, 
if this were not the case, this would hint towards severe 
truncation artifacts, and such problems could then be 
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FIGURE 3. The quenched quark mass function (top panel) 
and wave function renormalization (bottom panel) for the stan- 
dard gauge choice and the scaling gauge used in figure 1 from 
Dyson-Schwinger equations [16]. 

used to improve the truncations. This approach can be 
supplemented by more general considerations, as in QED 
[19]. However, the particulars of the non-Abelian case 
unfortunately makes this somewhat more intricate [20]. 

CONCLUDING REMARKS 

By using gauge-fixed correlation functions as a com- 
mon basis it is possible to combine different genuine 
non-perturbative methods in such a way as to extend 
their reach to problems none of them could reliably 
access alone. The basis for this combination has been 
constructed during the recent years. With these tech- 
nologies in place, this approach can now be applied 
to increasingly complicated problems. Prominent exam- 
ples are light, dynamical fermions [11] and the QCD 
phase diagram at both finite temperature and density 
[12, 13, 14, 15]. Furthermore, first investigations have 
started to also extend their reach to physics question be- 



yond QCD, in particular the electroweak sector in and 
beyond the standard model [6, 18, 21, 22]. Thus, this ap- 
proach to treat gauge theories beyond perturbation theory 
is one of the possible avenues to study strongly interact- 
ing systems at all scales. 

ACKNOWLEDGMENTS 

This work has been supported by the FWF under grant 
number M1099-N16. 



REFERENCES 

1. R. Alkofer, and L. von Smekal, Phys. Rept. 353, 281 
(2001), hep-ph/0007355. 

2. C. S. Fischer, A. Maas, and J. M. Pawlowski, Annals 
Phys. 324, 2408-2437 (2009), 0810.1987. 

3. A. Maas (2010), arXiv: 1010 . 5718. 

4. A. Maas, Phys. Lett. B689, 107-111 (2010), 
0907.5185. 

5. I. L. Bogolubsky, E. M. Ilgenfritz, M. Miiller-Preussker, 
and A. Sternbeck, Phys. Lett. B676, 69-73 (2009), large 
lattices, 0901.073 6. 

6. A. Maas (unpublished). 

7. A. Cucchieri, A. Maas, and T. Mendes, Phys. Rev. D77, 
094510 (2008), 0803.1798. 

8. W. Schleifenbaum, A. Maas, J. Wambach, and R. Alkofer, 
Phys. Rev. D72, 014017 (2005), hep-ph/0411052. 

9. A. Cucchieri, A. Maas, and T. Mendes, Phys. Rev. D74, 
014503 (2006), hep-lat/0605011. 

10. R. Alkofer, C. S. Fischer, F. J. Llanes-Estrada, and 
K. Schwenzer, Annals Phys. 324, 106-172 (2009), 
0804 .3042. 

11. C. S. Fischer, J. Phys. G32, R253-R291 (2006), 
hep-ph/0605173. 

12. D. Nickel, J. Wambach, and R. Alkofer, Phys. Rev. D73, 
114028 (2006), hep-ph/0603163. 

13. C. S. Fischer, A. Maas, and J. A. Miiller, Eur. Phys. J. 
C68, 165-181 (2010), 1003.1960. 

14. J. Braun, H. Gies, and J. M. Pawlowski, Phys. Lett. B684, 
262-267 (2010), 0708.2413. 

15. J. Braun, L. M. Haas, F. Marhauser, and J. M. Pawlowski 
(2009), 0908 .0008. 

16. M. Blank, A. Krassnigg, and A. Maas (2010), 
1007.3901. 

17. C. S. Fischer, P. Watson, and W. Cassing, Phys. Rev. D72, 
094025 (2005), hep-ph/0509213. 

18. A. Aguilar, and J. Papavassiliou (2010), arXiv: 
1010.5815. 

19. A. Kizilersu, and M. R. Pennington, Phys. Rev. D79, 
125020 (2009), 0904.3483. 

20. L. von Smekal, A. Hauck, and R. Alkofer, Ann. Phys. 
267, 1 (1998), hep-ph/9707327. 

21. A. Maas (2010), 10 07 . 072 9. 

22. F. Sannino, Acta Phys. Polon. B40, 3533-3743 (2009), 
0911.0931. 



